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We show that the vibrational state tailoring method developed for molecular systems can be 
applied for cold atoms in optical lattices. The original method is based on a three-level model 
interacting with two strong laser pulses in a counterintuitive sequence [M. Rodriguez et al., Phys. 
Rev. A 62, 053413 (2000)]. Here we outline the conditions for achieving similar dynamics with single 
time-dependent potential surfaces. It is shown that guided switching between diabatic and adiabatic 
evolution has an essential role in this system. We also show that efficient and precise tailoring of 
motional states in optical lattices can be achieved, for instance, simply by superimposing two lattices 
and moving them with respect to each other. 
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I. INTRODUCTION 

Complete control of atomic center-of-mass motion is 
one of the main principles of physics with cold atoms .1, 
0] and ions 0, 3 ■ Cooling of atomic samples allows their 
trapping with far-off resonant optical traps or with 
spatially inhomogeneous magnetic fields [g- This has 
made it possible, for example, to explore the quantum 
statistics of bosonic and fermionic atoms |B,l3|i to develop 
ultrahigh precision spectroscopy and atomic clocks 13 , 
and to produce cold molecules by photoassociation l9j. 
One specific example is an optical lattice, where the in- 
ternal state structure of atoms and the polarization states 
of light combine to produce periodic and controllable po- 
tentials for atoms J^]. By sufhcient cooling the atoms 
can be localized into the lattice sites. The situation is 
reminiscent of solid state systems, and interesting effects 
such as formation of a Mott insulator [Til IT^ or the sur- 
face of a Fermi sea ^3,^^ have been observed. This has 
motivated ideas of building quantum simulators for solid 
state phenomena, and even quantum computers 0, . 
For such purposes one needs a wide-ranging and versatile 
toolbox for the quantum control of atomic motion. 

In optical lattices atoms can either move in the Bloch 
bands of the periodic structure, or become localized at 
lattice sites (not forgetting the interesting intermediate 
region where the motion or "hopping" takes place by 
quantum tunnelling between the lattice sites). Here we 
consider the case of rather deep localization, where the 
eigenstates of the atomic center-of-mass states are quite 
discrete in energy and become the vibrational states of a 
single lattice site. Our aim is to develop a tool for mov- 
ing atoms efficiently and selectively to these vibrational 
states if they are initially in the vibrational ground state. 
Especially, we want to consider the case where the state 
change is associated with a move from one lattice site to 
another. 

We approach the problem with ideas originally pro- 
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posed for the quantum control of vibrational states in 
dimer molecules ^3 Note that for atoms the 

vibrational states are center-of-mass states, whereas in 
molecules they correspond to the internal motion, i.e., 
relative motion of the nuclei. An efficient method to move 
molecules between vibrational states is STIRAP (stimu- 
lated raman adiabatic passage) [2^], which is based on 
the existence of a noninteracting eigenstate in a three- 
level system and a counterintuitive pulse sequence. The 
advantage of the method is its high efficiency (popula- 
tion transfer is nearly complete), robustness in relation 
to the parameters, and the elimination of the intermedi- 
ate state. For molecules as well as for atoms in a lattice 
the intermediate state would usually be an electronically 
excited one, which normally decays rapidly and incoher- 
ently. 

The standard STIRAP approach, however, is strongly 
dependent on overlap integrals between the vibrational 
states, which in molecules means favorable Franck- 
Condon factors. This limits the achievable "stretch- 
ing" of the molecules as the process needs to adhere to 
the Franck-Condon principle. In addition, any imple- 
mentation for optical lattices would require additional 
pulsed laser fields. An interesting and viable alterna- 
tive in molecules is offered by very strong pulses, which 
actually affect the vibrational state structure |2j|. As 
demonstrated in Ref . (l^ , this situation can be modelled 
with LIP's, i.e., hght-induced potentials, which are the 
time-dependent potential surfaces associated with the in- 
stantaneous eigenstates of the combined system of three 
electronic states and two pulses. In fact, the process is de- 
scribed by adiabatic evolution along a single LIP, which 
efficiently takes the molecule from one electronic state 
to another. As the equilibrium positions associated with 
each electronic state may be quite different, the process 
leads to an apparent breakdown of the Franck-Condon 
principle. Later it was shown that with a proper set- 
ting of parameters one can decide whether the vibrational 
state is changed during the process as well . 

Here we show that the same method is applicable for 
atoms in optical lattices or in double- well structures. 
First we demonstrate that under suitable conditions we 



2 



can indeed reduce the system of three states and two 
pulses into a single time-dependent potential. As ex- 
plained in Sec. this is not so straightforward as one 
could assume. The fact that the process is not completely 
adiabatic (as noted in Ref. 18j) plays a role here, as well 
as in tuning the parameters to achieve the population 
of a specific vibrational state. In Sec. IIIII we provide 
some numerical results that demonstrate the practicality 
of the approach, and in Sec. lIIIEl we address specifically 
the two-dimensional case (the original molecular studies 
were performed in one dimension). A special case of two 
optical lattices that can be moved with respect to each 
other is treated in Sec. lIVI Finally, we conclude our pre- 
sentation in Sec. Ivl with a discussion on applications of 
the method. 



II. THREE-STATE AND SINGLE-STATE 
DESCRIPTIONS 

In case of molecular dimers, the system is originally 
described in terms of a three-component state vector con- 
sisting of electronic state wave functions Next we 
present a highly simplified picture of a dimer system. 
The advantage is that the description of the system can 
be generalized beyond the dimer concept to any suit- 
able three-state case, and then reduced to a single-state 
model. Both cases, three state and single state, are in 
principle applicable to suitable optical lattice systems, 
and furthermore, the latter model can be also used to de- 
scribe situations which are genuinely single-state cases, 
i.e., there is no underlying three-state structure present. 

Let us denote the state vector components ipii'^j t)^i = 
1,0,2 (following Ref. Consequently, the time evo- 

lution of such a state vector ^(x, t) is given by the scaled 
time-dependent Schrodinger equation 

»^*(x, t) - -V2*(x, t) + V(x, t)*(x, t), (1) 

where x and t are the scaled position and time. Because 
of the used scaling, every quantity is expressed with a di- 
mensionless number, and especially h = 1 and m = 1/2. 
The associated electronic potentials and state couplings 
form the matrix 

/i/i(x) ni{t) \ 
v(x,<) = n^it) T/o(x) n2{t) , (2) 

\ n2{t) F2(X) J 

where 

F,(x) = ic^f(x-xO' + A„ z = l,0,2, (3) 
n,{t) = n'^exp[~{t-tjf/Tf], J = 1,2. (4) 

In other words, we assume either harmonic or flat poten- 
tials and Gaussian pulses. Also, we have above applied 
the rotating wave approximation (RWA), which allows 
us to shift potentials appropriately. Thus the potentials 




Figure 1; (a) The initial system consists of three individual 
potential surfaces coupled to each other via two pulsed lasers 
with frequencies i/i. (b) In the rotating wave approximation 
(RWA) the system is characterized by Rabi frequencies fJi 
and detunings A;, (c) Finally, the light-induced potentials 
(LIP's) are derived from the RWA potential matrix, Eq. 
as its time and position dependent eigenvalues, shown here 
by dark lines. 



Vi(x) are parameterized by trapping frequencies uji and 
detunings , and the coupling terms flj (t) include max- 
imum Rabi frequencies f2j and time scales for the pulses 
Tj. The role of parameters and the effect of RWA is 
demonstrated in Fig. ^ 

We assume that Ai = and loq = 0. This assumption 
is justified because the effect of Ai is merely a universal 
shift in the energy levels and the shape of the potential 
Vb does not have any significant effect on the results 0| . 
We are interested in adiabatic time evolution with the 
exception that time evolution may include sudden rapid 
diabatic behavior. The potential matrix Q can always 
be diagonalized and we obtain three real eigenvalues as a 
result. Note that these eigenvalues depend both on time 
and position; they can be considered as new potentials in 
the eigenstate basis. We label these light-induced poten- 
tials (LIP's) with Vi(x, t). The eigenstate basis turns out 
to be equivalent to the adiabatic basis and the original 
three-state basis is called the diabatic basis. The num- 
bering of the LIP's goes as 1, 2, and 3, with increasing 
energy. 

In our model the couplings in the diabatic basis 
are initially (for t ^ ii,i2) and finally (for t ^ ^1,^2) 
practically but not quite zero. Then the diabatic and 
adiabatic potentials are very much alike, except that any 
degeneracy is a level crossing for the first case, and an 
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Figure 2: The time evolution of the light-induced potential 
surface Vi(x,f). For clarity, the energy values have been cut 
from above. Here we have Ao > and thus the flat poten- 
tial surface is located above the centers of the two harmonic 
ones, as in Fig. ^c). Due to the counterintuitive pulse order 
{ti > fa), the right-hand-side well drops down first, but the 
initially large position separation keeps the system in the di- 
abatic state 1 long enough. The grey dashed lines mark xi (t) 
and X2 (t) . 



avoided crossing for the second case, see Fig. ^c). If 
we are initially on the diabatic state 1, in the molecular 
model this means that adiabatic evolution will keep us in 
the appropriate electronic eigenstate, which at the final 
point turns out to have evolved into the diabatic state 2. 
For a system with three states this works very well, but 
if we merely take the eigenstate potential Vi (x, t) and 
solve the appropriate single-state Schrodinger equation, 
we find that initially our system corresponds to a double- 
well potential [combination of Vi(x) and V2(x); here xi 
and X2 are assumed to be clearly different] . As the right- 
hand well begins to move down in energy initially due 
to the strengthening coupling pulse 512, the separating 
energy barrier is too weak to hold the initial state lo- 
calized to the potential well Vi (x) , see again Fig. ^c) . 
In the three-state model this problem does not arise, be- 
cause initially the system is held by the diabatic potential 
Fi(x). 

In other words, the apparent adiabatic evolution along 
a single LIP is possible in the model only because in 
reality the evolution is diabatic initially and finally, and 
adiabaticity sets in only in the middle of the process. In 
order to work with truly single-state models, we need to 
modify the molecular approach. An obvious solution is 
to make Xi and X2 time dependent as well, by setting 
their initial and final values sufficiently large, so that the 
double- well potential barrier is strong when the adiabatic 
description would otherwise fail, see Fig. |21 



III. NUMERICAL RESULTS FOR THE 
SINGLE-STATE MODEL 

A. Description of numerical methods 

We consider the single-state wave function as a wave 
packet of plane waves and evolve it in a discretized space 
numerically with the methods of wave packet dynam- 
ics jl^l- We solve the time-dependent Schrodinger equa- 
tion using the split-operator and fast Fourier transfor- 
mation (FFT) methods. In addition, for a better under- 
standing of dynamics we follow Ref. and solve also 
the time-independent Schrodinger equation for fixed mo- 
ments of time using the Numerov algorithm (23i | . We 
simplify our notation by defining 



V'(x,t) 
F(x,t) 



iVi{^,t), ifA2>0, 

|F2(x,i), ifA2<0. 



(5) 
(6) 



Obviously, it would be seemingly practical to set the 
three-state model completely aside instead of using it as 
a tool for obtaining the single-state LIP, but with this ap- 
proach we can use the knowledge obtained with previous 
studies of the molecular system. 

To model the time dependence of Xi and X2, we rede- 
fine our diabatic potentials Vi and V2 as 



VS-(x,t) 
x,(t) 



[x-x,(t)p + A 



'3 ' 

-r Pj sech[7j(t - ts)], 
[ 0, aj + f3j sech[7j {t~ts)]), 



(ID) 
(2D) 



(7) 
(8) 



where ts is the moment when the individual potential 
wells are closest to each other (Fig.|21). For simplicity we 
use in our calculations equal maximum Rabi frequencies 
ill — = ^ ^■ii'i equal time scales of the pulses Ti — 
T2 = T. Pulses are applied in a counterintuitive order, 
that is ti > t2- We choose the time ts so that ts = 
^ (i 1 +i2 ) • Since the individual potential wells are initially 
apart from each other, we can consider them effectively as 
two separate harmonic potentials, and as a consequence, 
we select our initial state '(/'(x, to) to be case specifically 
the ground or excited state of the harmonic potential Vi . 



B. One-dimensional case 

In our scaled units time t evolves from values to 25 
and position ranges from —3 to 3. We choose the mo- 
ments for the pulse maxima as ti = 15 and ^2 = 10. 
The pulse time scales T are set to T = 5. The trap- 
ping frequencies of the harmonic potentials Vi and V2 
arc LUi = L02 = 20. As for the shifting parameters, we set 
ai = —0-2 = —1.5, /3i = — /32 = 1.0, and 71 = 72 = 0.30 
[cf. the trajectories of Xj(t) in Fig. [2]. As stated earlier, 
we have set = and Ai = 0. 
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Figure 3: (Color online) (a) Time evolution for a process that 
drives population from the state corresponding the left ground 
state into that of the right one (cf. corresponding potential 
surface in Fig-HJ. Here Ao = 850, A2 = 0, and fl = 274. (b) 
An alternative process where the excited vibrational states is 
reached. Here Ao = 850, A2 = -20, and Q = 162. 

The above settings leave as control parameters Aq , A2 , 
and fl. If we set Aq = 850, A2 = 0, and = 274, we 
obtain the basic case of a transfer from the ground vi- 
brational state of potential Vi to the ground vibrational 
state of potential V2 in the three-state model. The evolu- 
tion is plotted in Fig.|3Ja). The change from one ground 
state to another one is almost complete. 

As an alternative selection we can take Aq = 850, 
A2 = -20, and n = 162. Figure E^b) shows that the 
final state becomes the first excited state of potential V2 , 
as expected for the analogous three-state case. 

By further changing parameter values it is possible to 
achieve transfer from the ground state of Vi to various 
excited states of V2, as shown in Fig.2](upper panels). As 
in the three-state case, it is also possible to return com- 
pletely to the state Vi, but with a different vibrational 
state. The initial state can also be a vibrationally excited 
one for Vi. One can also make the transfer incomplete, 
which produces well-to-well superposition states. Note 
also, by comparing Fig.|3Ia) with Fig.^a), and Fig.|3{b) 
with Fig.^Jd), that the number of peaks at the midpoint 
of evolution can be different even if the initial and final 
states are the same. We return to this issue in Sec. 1111 151 

C. The role of parameters 

The number of parameters in the full three-state model 
is very large: maximum Rabi frequencies Sl^, detun- 
ings Ai , trapping frequencies Ui , time scales of pulses Tj , 
and pulse delay At — ti — t2. Furthermore, one can con- 
sider asymmetric pulse shapes, chirped pulses or some 
different potential geometry. Consequently, the param- 
eter space is too large for a complete mapping, and we 
limit our discussion to the effect of Aq, A2, and il within 



the previous parameter settings. With these three pa- 
rameters, we are already able to target arbitrary vibra- 
tional states of both Vi and V2 • 

The sign of Aq is essential. Choosing Aq < 0, i.e., 
locating Vo(x) below the other two curves, will result in 
the very robust APLIP process in which the wave 
packet is smoothly transferred from one ground state to 
another, but no other tailoring is possible. The opposite 
selection leads to much richer dynamics and excitations 
become possible as well but we pay for it with the 
dynamics becoming far more sensitive to the parameters. 
In short, as long as |Ao| is large enough, its sign is a suffi- 
cient switch between the two cases. From an alternative 
point of view, changing the sign of Aq is also equivalent 
to swapping the pulses fli <-*■ ri2- 

The role of Q is intuitively clear: the more powerful are 
the couplings, the more are the adiabatic potentials mod- 
ulated. On the other hand, and |Ao| counteract each 
other, such that the larger |Ao|, the larger values of are 
needed to modify the potentials by the required amount. 
Of course, they do not exactly cancel each other, since 
lower values of |Ao| tend to allow the wells to broaden, 
too. Nevertheless, the cancellation is efficient from the 
point of view of regions close to the minima of the wells, 
where the dynamics is also concentrated. Therefore, the 
counterplay between the detuning Aq and the coupling 
strength explains, e.g., the plotted stripe structure of 
Fig. 3 in Ref . |l3 • We have verified the existence of these 
stripes also in our single-state model. 

The detuning A2 between potentials Vi and V2 has an 
obvious function. If the potentials Vj are to be consid- 
ered one by one, they each have their own (and in this 
case the same) frequency uj. Because of our scaling, the 
corresponding energy difference between the harmonic 
eigenstates is AE = w in the scaled units. Starting with 
A2 = 0, we can find such combinations of parameters 
that will lead to the transfer of the population from one 
ground state to another as in Fig. OJa). The basic idea 
is that lowering V2 by A2 — — AE will drive the popula- 
tion now into the first excited state with approximately 
the same parameters [cf. Fig. Elb)]. This generalizes to 
higher excited states as well: A2 = — n x AE, n e N, 
will lead to the nth excited state of potential 2. Inter- 
estingly, by setting A2 = — (n-|- 1/2) x AE and keeping 

< 100 we always get the initial state, i.e., lowest vibra- 
tional state of potential 1 back, but by increasing O we 
can reach the higher vibrational states of potential 1. So 
the tuning process is clearly sensitive to the value of fl, 
as shown in Fig. |31 

In the mapping of Fig. we can see that the final out- 
come forms regular patterns. Especially, clear system- 
aticity arises. On the one hand, higher vibrational states 
of Vi are gained by increasing while A2 acts only in a 
periodic way. On the other hand, the higher vibrational 
states of V2 are achieved by decreasing A2 while fl has 
no significant role. 
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Figure 4; (Color online) Examples of time evolution (upper panels) and corresponding energy spectra of six lowest LIP 
eigenstates (lower panels) for different parameter values. We have chosen Aq = 850 in every frame while A2 and are varied, 
such that (a) A2 = -5, Q = 100, (b) A2 = -5, n = 140 ( c) A2 = -5, ^ = 200, (d) A2 = -15, Q = 100, and (e) A2 = -35, 
— 175. Other parameters are as mentioned in Sec. Ilil 51 In the energy spectra, the populated state is marked with arrows 
and the encountered diabatic jumps with circles. 
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Figure 5; (Color online) The final states as a function of Q 
and A2. In the labels _R and L signify whether the wave 
packet resides at the end in the right-hand or left-hand well, 
respectively, and the number indicates the vibrational quan- 
tum number of the final state. 



D. Interpretation in vibrational basis 

In order to understand the role of the parameter 
choices and the dynamics in general, we have followed 
Rcf. 18] and analyzed the situation by solving the in- 
stantaneous vibrational states of the relevant LIP (we call 
them LIP eigenstates). If the two wells are sufficiently 
separated, the lowest few LIP eigenstates are very much 
similar to those of two single harmonic potentials Vi and 
V2 , except for the narrow degenerated cases, when the en- 
ergies of two LIP eigenstates coincide. As the energy of a 
LIP eigenstate gets comparable to the height of the sep- 



arating potential barrier, it starts to have non-negligible 
component in both wells, as seen in Fig. 

During the pulses the double-well potential gets mod- 
ulated in such a manner, that the potential barrier be- 
comes lower or even disappears. Accordingly, the set of 
corresponding LIP eigenstates changes as well. As the 
height of the barrier gradually decreases, the minima of 
the two wells move with respect to each other in energy. 
With large enough values of the changes are consid- 
erable, if compared to the energy separations AE. Con- 
sequently, in terms of the LIP eigenstates, the states ap- 
proximately corresponding to the single-well eigenstates 
appear to cross each other. This phenomenon can be 
easily seen in the energy spectra (cf. Fig.0J). Of course, 
there is always a small energy gap between the energies, 
i.e., the crossings are in fact avoided. 

The sharp crossings are an indication of diabatic evo- 
lution. Let us define crossing time scale 

r,",„,, = min {<' > 0| |(0„(t)|^„+i(t + i'))l' - l}, (9) 

where (f>n{t) is the nth LIP vibrational eigenstate. With 
the given choice of parameters, the transformations from 
one eigenstate to another happen in time scale of t^qss ^ 
10~^, which is less than the vibrational time scale of the 
potential Tvib ~ 10~^. Accordingly, the wave packet can- 
not succeed in following these changes in the vibrational 
basis, and therefore a diabatic jump between the LIP 
eigenstates occurs. 

Between the sharp diabatic crossings, the system fol- 
lows a particular LIP eigenstate adiabatically. Adiabatic 
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Figure 6: (Color online) The evolution of the double-well potential surface and the corresponding LIP eigenstates. This series 
corresponds the choice of aprameters in Fig. Bid)- Each frame represents a time slice, such that (a) t = 0, (b) t = 7.5, (c) 
t = 12.5, (d) t = 18, and (e) t = 25. For clarity, the amplitudes have been rescaled and the populated state is plotted with a 
thicker line. A rapid diabatic jump has occured at t ~ 7. 



following occurs when the energy gap between the two 
approaching LIP eigenstates is large. Consequently, in 
case there are many possible energy level crossings, the 
time scale of potential modulations Tpot can be used as a 
switch to determine which crossings are passed diabati- 
cally and which adiabatically. By making a crossing into 
a mixture of diabatic and adiabatic, one can reach super- 
positions, although the precise result is likely to be very 
sensitive on parameter values. 

It is worthwhile to emphasize, that the seemingly rapid 
and radical changes in the appcarcnce of the wave packet 
typically around t ^ 10 and t « 15 do not have anything 
to do with diabatic crossings. In Fig. we have plotted 
the LIP eigenstates corresponding the choice of parame- 
ters in Fig. ^d). By comparing these figures, it is obvi- 
ous, that the diabatic jump occurs already at i « 7, when 
nothing seems to happen in the time evolution. Interest- 
ingly, the subsequent changes in the course of evolution 
are in fact due to the evolution of the corresponding LIP 
eigenstate itself. Therefore, all the visible transforma- 
tions are actually a manifestation of adiabatic following. 

Returning to the discussion of Fig. \^ in the preced- 
ing section, we are now ready to give explanation to the 
pattern. Let us denote the number of sharp LIP energy 
level crossings by (n,m), where n and m, respectively, 
count the crossings before and after the pulse maxima. 
For example, the situation in Fig.^fa) corresponds to the 
notation (0, 1). 

The edges of the zones arise for two reasons. First, 
those with a positive slope occur because a new LIP en- 
ergy level crossing is being formed prior to the pulses, i.e., 
(n, m) ^ (n -|- 1, m) at the edge. Second, those with a 
negative slope are associated with a creation of a crossing 
after the pulses, i.e., (n, m) (n, m + 1) at the edge. Yet 
another rule for crossings can be formulated, but it does 
not correspond to any visible edge in the mapping. The 
third rule is, that when passing the resonance detuning 
A2 = (mod AE) from above, one crossing travels via 
infinity from left to right, i.e., {n,m) {n — l,m + 1). 
Illustration of the rules is given in Fig. [7| (cf. spectra in 
Fig. 21 and the mapping in Fig.jSJ. 



(a)(n,m)^ / (b) ^{n,m + l) 

' (n + 1, m) (n, m)' 





(c) {n, rn) 
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Figure 7: The three possible types of changes in the numbers 
of sharp LIP energy level crossings before and after pulse max- 
ima (n,m), respectively, at the edges of the zones in Fig. |3 
(a) Positive and (b) negative slopes affect different numbers, 
(c) An additional rule associated with the resonance energies 
A2 = (mod AE) (dashed line; not a visible edge). 



With the three rules, we can see that |n — to| < 1. If 
n = TO, then the nth vibrational state of potential Vi is 
achieved. If n = to ± 1, then a state in V2 is reached, 
such that each increment of —/S.E = —20 in A2 intro- 
duces one step higher excitations. Advancing towards 
greater values of Q. brings along more diabatic jumps. 
Naturally, the stability suffers as the complexity of the 
process increases, which can be seen as a degeneration of 
the pattern with large values of fi. 



E. Two-dimensional case 

So far we have studied only the one-dimensional case. 
In two dimensions, we set the potentials symmetric for 
simplicity and apply a time-dependent shift to the poten- 
tials Vi and V2 . The shifting is the same as previously in 
the one-dimensional case, but it is conducted only in one 
direction, as defined in Eq. ||SJ|. 

There is no coupling between the Cartesian coordinates 
X and y in the potentials Vi and V2 ■ Accordingly, Eq. ^ 
is separable in x and y. Thus the two-dimensional system 
is in practice split into one-dimensional slices which are 
labelled by a; = constant. The only difference between 
two X = constant slices is a shift in energy. Therefore, 
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Figure 8: (Color online) The double- well time evolution in 
two dimensions. The frames show the stage of evolution at 
(a) t = 0, (b) t — 10, and (c) t = 20. The parameters are 
the same as for the corresponding one-dimensional case in 
Fig-lib). 



we can expect the system to behave in a similar way as a 
one-dimensional system. Our numerical results confirm 
this, see Fig. |S| 

We have used the same parameter values in our calcu- 
lations as in the corresponding one-dimensional case. It 
should be emphasized, though, that the two-dimensional 
system is much more sensitive to specific parameter val- 
ues than the one-dimensional system. Another difi'erence 
is the fact that we can apply the time-dependent shift 
to the potentials Vi and V2 only in one direction [23 |. 
Thus the possible excited state outcomes of the processes 
have excited vibrational states in that particular direc- 
tion only. For optical lattice configurations that opens 
interesting possibilities. For example, one can shift the 
potentials in y direction in the first process and then take 
the final state of that as the initial state of another pro- 
cess where the roles of y and x are reversed. 



IV. APPLICATION TO A MODULATED 
LATTICE 

A. The basic idea 

The preceding model discussed so far originates from a 
three-state model and therefore inherits a somewhat un- 
intuitive parametrization by pulse strengths and detun- 
ings. The actual geometry of the double-well potential is 
given as a solution of a third degree polynomial result- 
ing from an eigenvalue calculation of the matrix given 
in Eq. (0). The solution cannot be simplified into any 
illustrative and compact form in terms of the control pa- 
rameters, which is also why it is not given explicitly in 
the text, and mainly qualitative characterizations could 
be given in Sec. IIII CI Nevertheless, the previously de- 
rived one-state model is adequate for presenting the gen- 
eral behavior of wave packet dynamics in time-dependent 
double wells. Therefore, the results do not depend on the 



Figure 9: (Color online) The sinusoidal potentials V{x) and 
W{x,t) used in the lattice modulation model. Superposition 
of these two results in a lattice of modulated double-well po- 
tentials U{x, t),asW moves with respect to V. The strengths 
of the lattice potentials are chosen such that Vb = 2Wo = 30. 
The dot marks the adjustable lattice offset parameter xq. The 
frames represent the situation with (a) xq = and (b) xq — 6. 



actual formulation of the underlying potential. 

We now move on to seek alternative realizations for 
the double-well potential. The main focus is to make the 
control of the potential as self-explanatory as possible. 
Another point of interest is to introduce a method which 
could be also experimentally accessible. 

The previous dynamics included typically two sets of 
diabatic crossings: one before the pulse maxima and an- 
other one after them. We now propose a possible scheme 
which will include only one set of energy level crossings. 
With the interest in optical lattices, we formulate the 
problem in terms of sinusoidal potentials. The modula- 
tion of the potential is done by another sinusoidal poten- 
tial (but with a doubled frequency), which can be moved 
with respect to the other one (cf. Fig.^. Consequently, 
each site of the potential lattice will now become a time- 
dependent do able- well potential. 

Let the wavelengths of the potentials be Xi in scaled 
units. Then the respective potentials are 

V{x) = Vasi-a^{2TTx/Xv), (10a) 
W{x,t) = Wosin^{2TT[x-xo{t)]/\w}, (10b) 
U{x,t) = V{x)+W{x,t), (10c) 

where Ay = 2X\y and xo{t) describes the offset of W with 
respect to V. The initial position is chosen such that 
the minima of both V and W coincide, i.e., a;o(0) = 0. 
Furthermore, approximately Vq « 2Wo is needed. 

The scaling is chosen as previously with an added as- 
sumption that Ay = 2Xw = 40. Accordingly, the recoil 
energies i?rccoii — 47r^ft^/(2mAf ) of the lattice fields be- 
7rV400 and E^^^.^ = TrVlOO in the scaled 



come EY^^^,^ = 
units. 

For simplicity, let us consider only one lattice site lo- 
cated initially at a; = 0. At first, the wave packet is at rest 
at the bottom of the well. Thereafter, the W potential 
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Figure 10: (Color online) The energy spectrum during lattice 
modulation as a function of the lattice offset parameter xq. 
Here Vo = 30 ^ 1220 EY^^^^ and W^o = 15 152 E^^^^ 
are used. The initial state is marked with an arrow. Using 
a sweep rate as in Fig. 1111 a superposition of the first two 
excited states (grey arrows) is achieved by passing the first 
crossing diabatically (black circle) and the second partially 
diabatically (grey circle). 



sweeps one-half of its wavelength to the right. Mean- 
while, the wave packet gets driven uphill by the moving 
potential barrier. However, since the barrier effectively 
becomes lower after having passed the a; = point, the 
wave packet will drop down to a vibrational state of a 
newly created well around x = 0. 

The point is that from the perspective of the initial 
ground state there will be only one set of energy level 
crossings during the period, when the wave packet is 
driven upwards in a local minimum, while a new global 
minimum is deepened on the other side of the barrier. 
The energy gaps in these crossings increase with the cor- 
responding vibrational state number (cf. Fig. ll()|l . There- 
fore, the sweep rate Tpot will determine which crossings 
are passed diabatically, and which are followed adiabati- 
cally. 



B. Numerical results 

Sinusoidal potentials are relatively flat and round if 
compared to the LIP's described in Sec.^ Therefore, the 
corresponding double- well structure holds a set of eigen- 
states, which do not represent eigenstates of individual 
wells very well, but are spread over both wells instead. 
Consequently, in case of the higher vibrational states the 
energy level crossings are not as sharp as previously. Ac- 
cordingly, the population is after the sweep likely to be 
spread over a distribution of vibrational states, as the 
population leaks gradually at each energy level crossing. 
Alternatively, a high control over the sweep rate Tpot{t) 
is needed in order to slow down at the instant when adi- 
abatic following during a crossing is wanted, and accel- 
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Figure 11: (Color online) Time evolution during the lattice 
modulation. The second plot is a high-resolution view of the 
transition moment, and illustrates the speed of the process 
compared to the smoother behavior in the harmonic wells. 
The lattice parameters are chosen as in Fig. llUl and a constant 
sweep rate is applied. 



erate again thereafter in order to cross the other possible 
gaps diabatically. 

Figure shows an example of the lattice evolution 
corresponding to Fig. ^] The lattice offset is evolved 
linearly, i.e., xo{t) — t/rpot, where the constant sweep 
rate is Tpot = 10**. A superposition of the first two ex- 
cited vibrational states is populated accordingly. The 
above-mentioned leaks of population to other vibrational 
states result in small oscillations of the final wave packet. 
The process of change is rather fast as expected, too. 
Nevertheless, this example shows clearly that the ideas 
developed for the three-state system with two pulses can 
be carried over to the optical lattices. 



V. DISCUSSION 

We have investigated a transition from a counterin- 
tuitively coupled three-state system into a single-state 
system described by a time-dependent double-well po- 
tential. Thereafter, we have shown how the populations 
of the vibrational states in each well can be tailored by 
modulations of this double-well potential. 

The focus of our study was in investigation of necessary 
diabatic crossings blended among otherwise adiabatic fol- 
lowing. The system holds basically three different time 
scales: (i) the intrinsic vibrational time scale Tvib induced 
by potential geometry, (ii) the general time scale of the 
modulation of the potential surface Tpot , and (iii) a set of 
time scales associated with the crossings of the LIP eigen- 
states 71 G N with their upper neighbors {t^q^j,}. First, 
Tvib <C Tpot is needed in order to have adiabatic following 
in general. Second, the tailoring of the state populations 
is done by adjusting the set {t^ossI with respect to Tvib 
such that certain energy gaps are crossed diabatically 
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while the others are avoided adiabatically. The values 
of T|!^oss are ultimately determined by the way the po- 
tential geometry is modified, but Tpot scales the overall 
time scale of the process. Therefore, the adjustment of 
crossing time scales can be made by tuning either the 
potential parameters or the general time scale Tpot • 

In this general model survey we have used scaled units 
rather than taking a specific physical system with fixed 
parameters. This was done intentionally in order to fo- 
cus on describing the model itself. In practice, we expect 
that the results can be applied in various systems. A 
recent example is an optical lattice setup where one gen- 
erates a periodic two-dimensional lattice of double-well 
structures, with controllable barrier heights and relative 
well depths '2^1 . 

Another interesting situation arises when we consider 
our model from the viewpoint of a moving atom in an 
atom chip waveguide |26j| . The position in our descrip- 
tion can be considered as the transverse direction of the 
waveguide, and the time evolution corresponds to the 



change in the potential seen by the atom during its mo- 
tion in the longitudinal direction. Our model can then 
describe a controlled transfer process between two waveg- 
uides, with the additional possibility of exciting tranverse 
modes. 

Our future investigations will consider the role of inter- 
actions in this process. The simplest system to study is 
a zero-temperature Bose-Einstein condensate. The inter- 
actions cause the LIP eigenstate structure to depend on 
the momentary local density and this back-action needs 
to be taken into account. Also, in two dimensions it will 
remove the separability of the coordinates. 
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